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Abstract 

Higgs pair production at the LHC from gluon fusion is small in the Standard Model, but can be 
enhanced in models where a resonant enhancement is allowed. We examine the effect of a resonant 
contribution from a second scalar arising in a model with a gauge singlet scalar field in addition 
to the usual SU{2) scalar doublet, with mass up to Mh ~ 600 GeV and discuss the interference 
effects in double Higgs production. The interference effects distort the double Higgs invariant mass 
distributions, and, depending on Mh, can enhance the total cross section by up to ~ 20% or 
decrease by ~ 30% for viable mixing parameters. We compute the NLO QCD corrections in the 
large mt limit. The corrections are large and can also significantly distort kinematic distributions 
near the resonance peak. 
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I. INTRODUCTION 


The experimental exploration of the Higgs sector of the Standard Model (SM) is one of 
the main goals of the current LHC run. Current data on Higgs properties are in reasonable 
agreement with the theoretical expectations, although there is still considerable room for 
new physics. An attractive extension of the SM is the Higgs portal scenario, in which the 
SM Higgs boson couples to a gauge singlet scalar, S, which in turn can communicate with 
a hidden sector. Models with an additional scalar singlet have also been used to generate a 

nn 

strong first order electroweak phase transition[lH^. 

In the Higgs singlet model, the SM Higgs doublet mixes with the new singlet, S, to form 
two physical scalar bosons: one, h, identified with the observed rrih = 125 GeV resonance 
and a second, H, with mass Mh- When Mh ^ 2m/i, large resonant enhancements are possi¬ 


ble in double Higgs production from gluon fusion, significantly enhancing the rate compared 
to the SM prediction. The singlet model has the advantage of depending on relatively few 
parameters, allowing for straightforward experimental study at the LHC in the analysis of 
Higgs couplings j? I, searches for heavy SM-like Higgs bosons 8-10| and direct searches for res¬ 


onant di-Higgs production |lll-ll4|. Higgs singlet models have also been extensively studied 


theoretically and additional limits derived from precision electroweak data, the interpreta¬ 


tion of LHC results, and restric 
perturbativity of the couplings 


ions from the requirements of perturbative unitarity and 
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3Q|. 


Double Higgs production from gluon fusion in the SM results from both triangle and 
box loop contributions, which interfere destructively, causing a suppression of the total rate 
from the naive estimate}^, 32|. This process has been studied at lowest order QCD (LO) in 
the singlet model, and regions of parameter space with enhanced rates determined. In this 
work, we consider precision predictions at NLO QCD for double Higgs production in the 
singlet model, including the hh invariant mass distribution. Since double Higgs production 
from gluon fusion first occurs at one-loop, the full NLO corrections involve two-loop virtual 
diagrams with massive internal particles. The calculation is considerably simplified by using 
an effective theory corresponding to the rrit ^ oo limit of the SM. In the SM, the corrections 
to the total rate have been known at NLO for some time in the effective theory 33(], which 


has also been matched onto t 
been calculated at NNLO 35 


le IN 


NLL threshold resummed result 


34 


36| and matched to the NNLL result 


37 


. lecently the rate has 
These corrections 
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typically increase the rate by a factor of about 2 — 2.3. The SM NLO QCD corrections to 
gg —)■ hh are also known in an effective held theory limit where the exact mass dependence 
is retained everywhere except in the virtual corrections 38j and alternatively in an expansion 


m 
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4n| . The unknown rrit dependence of the higher order QCD corrections induces 


an uncertainty of O{±10%) in the SM predictions. 

Higher order QCD corrections to new physics scenarios with resonant enhancements of 


3,0 


the double Higgs rates have been derived for the MSSM 
model 10, and also in an effective operator formalism with 


and the two 


no resonance 


diggs doublet 


43] . These cor¬ 


rections not only affect the total rate, but in some regions of parameter space distort the 
shape of the distributions. In this paper, we examine the approximations behind the QCD 
corrections in the context of the Higgs singlet model. We demonstrate that the corrections 
in the resonance region are signihcant and that the use of a constant K factor is a poor 
approximation in this regime. We also investigate the interference effects between the heavy 
scalar and SM-like contributions. These effects can be signihcant and should be included in 
searches for new heavy scalars. 

II. MODEL 


A. Recap 

We consider a simple extension of the SM containing the SM Higgs doublet, <I>, and an 
additional real gauge singlet scalar, S. After imposing a Z 2 symmetry under which S ^ —S, 
the most general scalar potential is [3, 171 

/In ^ JIa . 

( 1 ) 


V = - 11 ^ #t<I) - + jS‘. 


s without 


6|. How- 


Although not necessary for a strong hrst order electroweak phase transition, moc 
a Z 2 symmetry have been constructed in the context of electroweak baryogenesis 
ever, the additional complication is not necessary for our discussion of higher order correc¬ 
tions. After spontaneous symmetry breaking, in the unitary gauge we have = (0, (j)o)/ 
with (0o) = V = 246 GeV and S = {s + x)/\/2 with (S) = x/'/2. 

The mass eigenstate helds, h and H, are: 

h \ I cos 9 — sin 6^1 / (po — v 

H ] \ sin6' cos 9 / \ s 


( 2 ) 
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with physical masses, rrih and Mh, and — f < 6 * < |. 

The terms in the potential can be written in terms of the physical masses and mixing 
angle as, 

1 


2 2 \ I 2 

/i = V X + -X a2 


m 


= - {x‘^b 4 + ^^ 02 ) 


rni , Mf, 


^ ^ + 


mt 


2n2 


sin^ 9 


0-2 — 




mt 


vx 


(2 sin 9 cos 9) 




+ 


2 (m 


I-Ml) . 2 ., 

^ —— sm^ 9 . 


x^ x^ 

The requirement that the potential be bounded from below imposes. 


02 > —2^/b^ , 

We will also need the triple scalar couplings: 

L 


A, 64 > 0 . 




where 


Am — 


V 


6 


cos^ 9 — tan (3 sin^ 9 


m? f A 

A211 =-^ sin 26 '(cos 0 + sin^tau/S) ( 1 H -^ 

V \ 2miJ 


( 3 ) 

( 4 ) 

( 5 ) 

( 6 ) 

( 7 ) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 


and tan/? = -. A complete list of the scalar self-couplings can be found in the Appendix of 

Ref. Q- 

We assume that the lightest scalar, h, is the SM-like Higgs particle with irih = 125 GeV. 
The decay widths to SM particles, X, are then simply the SM values rescaled by the scalar 
mixing angle, 

T{h XXl) = cos^ 9T{h XX^)sM 

T{H XXl) = sin^ 9T{H XX^)sm (12) 

where r{H —)■ XX'^)sm is the SM partial width evaluated at mass Mh- The total widths 
are 


W = cos^^T 


SM 


Th = sin^ + 9{Mh - 2mH)T{H hh), 


( 13 ) 
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FIG. 1: Branching ratio for H —?■ hh. 


where is the SM total width evaluated at mass and 


hh) = 


^211 


Ami 


Ml,- 


(14) 


327iMh 

The branching ratio of H ^ hh is shown in Fig. [1] For small sin0, the branching ratio is 
relatively insensitive to tan/3 and is approximately BR(i7 —)■ hh) ~ 0.3 — 0.4. 

The model has 5 free parameters which we take to be: 


mh = 125 GeV, Mh, v = 246 GeV, tan /3, cos 6 . 


( 15 ) 


B. Limits 


The Z 2 symmetric Higgs singlet model is restricted by a number of experimental mea¬ 
surements. Fits to the h couplings assuming no branching ratio to invisible par ticles require 
I cos 6 \> 0.93 at 95% conhdence level[^. Precision electroweak quantities]^, in particu¬ 


lar the W boson mass45|, receive contributions which are sensitive to Mh and cos 6 *. For 
Mh ^ 400 GeV^ measurements of the W mass require | cos9 \> 0.96, with the limits signih- 


cantly weaker for smaller Mh 15 


3,Q- 


Heavy Higgs searches can also be interpreted as limits 


on cos 6 *. For Mh ^ 300 GeV, these limits are stronger than the limits from the W mass. 
Assuming no branching ratio, H —)■ hh, the direct search limits for heavy Higgs bosons can 
be interpreted as requiring | cos 6 \> 0.92 in this region. Requiring 64 to remain perturbative 
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FIG. 2: Feynman diagrams contributing to gg —)• hh in the singlet model. 


as it is scaled to high energy gives an upper limit on tan (3 w 


sinO = 0.1 and Mh = 200(500) GeV, tanf3 < 1.5(0.5) ld,l2Q[. With these considerations in 


rich depends on Mh and 9\ for 


mind, we will in general present results with cos 6* = 0.96, tan/3 = 0.5. 


III. DOUBLE HIGGS PRODUCTION 


A. LO Results 


Two Higgs production arises from the Feynman diagrams shown in Fig. [2l The result is 
sensitive to new colored objects with mass m (fermions or scalars) in the loops 46l-l5l| and 
also to the 3— Higgs self-couplings. The amplitude for t h{p')h{q') can be 

written as. 


= 

87rn2^ 


Sab{ Pr{p,q)F'i{s,t,u,m^) + P^’'{p,q,p')F2{s,t,u,m^ 


( 16 ) 


where Pi and P 2 are the orthogonal projections onto the spin-0 and spin-2 states respectively, 

PriP.Q) = 

p-q 

Pn?^ Q,p') = 9^'" + ^{+ {t- mDq^p'’' + {u - ml)p''p'>^ + spV'| , (17) 

spr f J 

s, f, and u are the partonic Mandelstam variables, 

s = {p + qf 


t = (p-p) 


/\2 


u = {p-q') 


/^2 


Pt = 


ut — 


(18) 


The functions Fi and F 2 are known analytically!^, 321 and the partonic cross section is 
given in terms of the form factors by. 


d^To ^ ( \ p + I F2{s,t,u,m^) 

dt \ 


2 \ 12 


(19) 
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where hr is the renormalization scale. (We have included the factor of | for identical 
particles in the hnal state). In the singlet model (as in the SM), the dominant contribution 
comes from top quark loops. The form factors can be written as, 


= —s 


cos6*Aiiin 


+ 


Fi(s, t, u, ml) = Ff *(s, t, M, ml) + t, u, ml) 

Fl^^{s,t,u,ml) 

Fi°^{s,t,u,ml) 

F2{s,t,u,ml) 


sin6*A2iin 


s - ml + imhTh s - Mfj + iMnT h 

s cos^ 6*Fn(s, t, u, ml) 

s cos^ 6*Gn(s, t, u, ml). 


F^{s,Tnl) 


( 20 ) 


In the limit —?> cx), 


F/\ —)■ 


Fn 


4 

3 


Gn —)■ 0 . 


( 21 ) 


The form factors Fa ,Fn, and Gn including the full kinematic dependences are found in 


Refs. 


31 


32|^. We denote the cross section found by including the exact mt dependence of 


the matrix elements, Eq. [201 by a'^Q, and the mt ^ oo limit, Eq. EH as ■ 

The LO hadronic cross section is. 


<6= l^^dT^d-S{s = TS), 


and the luminosity function is defined. 


dr 


ij F ^ 


/iF 


S is the square of the hadronic energy, tq = and /if is the factorization scale. 


( 22 ) 


(23) 


B. NLO Corrections 


The NLO corrections in the SM are known in the large mt limit 33|] and are trivially 
generalized to the singlet model. The gg initial state contains IR singularities which cancel 


^ The functions defined in Eq. ED] satisfy EAlEbjGn) —>■ 2FA(Eb,Gn)(Ref. 311) and sEA(sF:,sGn) 

i^A(Fn,Gn)(Ref. 0)- 
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when the real and virtual contributions are combined. The remaining collinear divergences 
in the gg, qg and qq initial states are absorbed into the NLO PDFs defined in the MS 
scheme with 5 light fiavors. The terms listed below are the finite contributions obtained 
after canceling the singularities. We write the NLO rate as, 


'^NLo(PP - ^LO + ^Tirt + ’ 




(24) 


where, 


cr""‘ = 

‘^virt 


= 

99 


(Jj / 


TT 


'TO 




dr 


71 


'TO 


dr 


'tq/t 


— { -^^99(^) log- 


Mf 


TS 


11 


—{i-zY + 6[i + z^ + {i-zy] 


/ log(l - z) 

V 1-^ 


= 

99 




TT 




»log 




rs(l — zy 


+ -z^-{i-zr 


= 

99 


OisifJ^R) 


TT 


f\ ^dC^^ dz 32, 


(25) 


We follow the philosophy of Ref. 33| and approximate the form factors in the virtual 
corrections by the exact dependent quantities and include the full mass dependence in 
d^Q in Eq. [25l The coefficient, C'™'‘, for the virtual corrections is 

= ^2^11^ 33-2n,; 

2 6 s 


85 

— cos^ 9Real 
9 




dtl Fi{s, t, u, ml) - ^{s - 2 ml)F 2 {s, t, u, mj) 


^ dt^ \ P + I F2{s,t,u,ml) 


(26) 


and 


,/i-Li 


(27) 


Pgg{z) and Pgq{z) are the DGLAP splitting functions. 


^ 99 (^) = 6 


1 \ 1 . s 

— ] +--2 + z[l-z] 


33 — 2n 


If 


6 


<5(1 - z), 


^ 4 i+(i-zr 
”< ■' 3 2 


( 28 ) 





















where nif = 5. The result in Eq. [25] has only approximate hnite rrit dependence since it has 


been adapted from the NLO calculation in the m* —)■ cx) limit 
We dehne an rrit dependent differential K factor from Eqs. 


33| 


and 


— ^^NLO 


daZ 


(29) 


dMhh / dMhh 

where Mhh is the invariant mass of the final state double Higgs system. In the following 
section, we will also show the numerical effects on the K factor of replacing the form factors 
and LO cross section by their rrit ^ oo limits, 


NLO 


j^mt^oo — 


dM, 


hh 


dM, 


hh 


(30) 


where. 


2 o Read ca{s) - cos"^ 9 

C-”-- = + T + - 

2 D so I CA(Sj — COS^ p 


and 


ca = 


cos 6'Aiiin 


+ 


sin 9X2uV 


s - ml + imhVh s - Mfr + iMuT 


H 


(31) 


(32) 


IV. RESULTS 


Our results are computed using CT12NLO PDFs|52| with a central scale choice hr = 
Pi? = /i = Mhh for the renormalization and factorization scales, and with mt = 173.34 GeV 


and mfo = 4.62 GeV. In the computation of T{H —)■ hh), we use the MS NNLO running 
mass for mh{MH) and we always assume that = 125 GeV is the lightest Higgs boson. 
Finally, the production cross section is computed including only the top quark loops, which 


are the 
HPAIR 


argest contribution. Our numerical results in the SM are checked using the program 


33[ |. The singlet model results from our private code were checked by incorporating 


the resonance from the singlet model in HPAIR and comparing the two results. 


SM Results 


The LO rate for gg —)■ hh in the SM is well known, as are the NLO and NNLO rates in 
the rui —oo limit. Ref. 37] hnds the NNLO matched to NNLL rate of 36.8 fb for pp hh 


at y/S = 13 TeV, p = Mhh, using MSTW2008 PDFs. The contributions to the differential 
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pp^hh (SM), Vs = 13 TeV pp^hh (SM), Vs = 13 TeV 



FIG. 3: Contributions to the SM rate for pp —)• hh at ^fS = 13 TeV in the —)• oo limit (LHS) 

and using the approximated mt dependence of Eqs. I25I26I (RHS). 

SM NLO cross section are shown in Fig. [3] in the oo limit (LHS) and in the 

dependent approximation of Eqs. [25] and [26] (RHS). The normalization and shapes of the 2 
approximations are quite different, but the K factors computed from the 2 approximations 
are almost identical. The contributions from real gluon emission, (jgg, and from the one-loop 
virtual diagrams, cXvirt, are of similar sizes, while the contributions from quark initial states 
are highly suppressed. In Fig. [4] we show the NLO result with approximate mt dependence 
as dehned in Eqs. [25l [26l and LO results for —)■ oo and including the dependence 

exactly. The lowest order result in the mt ^ oo limit overshoots the exact lowest order 
result at high M^h and fails to reproduce the peak structure, as is well known. Including 
the NLO corrections signihcantly increases the rate. (Calculating from the LHS of 

Fig. [3] and weighting by the exact mt dependent LO result gives a curve which is almost 
indistinguishable from the NLO curve of Fig. [T]) 

We show the renormalization/factorization scale variation of the SM LO and NLO rates 
in Fig. [5] when Mhh/‘2 < fi < 2Mhh- In this hgure, the LO rate includes all top mass 
dependence and the NLO rates are calculated using Eq. [25] The fractional scale dependence 
is signihcantly reduced at NLO. The scale variation of the differential SM factor dehned 
in Eq. [29] is shown in Fig. [6] At = 400 GeV^ the NLO scale uncertainty is ~ 11%, 
while at = 800 GeV it is ~ 15%. In the SM, the diherential K factor is only slightly 
dependent on Mhh and can be accurately approximated by a constant. 
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pp^hh (SM), VS=13 TeV 

CT12PDFS, r\=125 GeV 



FIG. 4; SM differential cross section for pp —)• hh at ^/S = 13 TeV. The NLO curve labelled cr™'* 
is obtained from Eqs. 1251261 


pp^hh (SM), VS=13 TeV 

CT12PDFS. m|=125 GeV 



FIG. 5: Scale dependence of the SM differential cross section for pp —)• hh at \/5 = 13 TeV. The 
NLO curves are obtained using Eqs. 1251261 

B. Singlet Model Results 

We begin by showing some lowest order results. The LO rate as a function of is 
shown in Figs. 0 and [H For the smaller values of tau/S and 6, the resonances become 
narrower, while for heavier the height of the resonance peak and the dip above the 
peak due to interference effects become smaller. The strength of the destructive interference 
is particularly strong for Mh = 200 GeV. Interference effects will be more thoroughly 
discussed in the next section, but we give an outline here. As in the SM, the box diagram 
dominates the h-resonance. Hence, the major contributions to interference are between 
the h- and iL-resonances and the box diagram. The h-resonance and box diagrams (SM 
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pp^hh (SM), VS=13 TeV 

M^j^/2 < ji < 2 CT12PDFS, m^=125 GeV 



M^JGeV) 


FIG. 6: Scale variation of the SM K™* factor, Eq. [29l for pp —)■ hh at ^/S = 13 TeV when the 
scale is varied < p < 2Mhh- 

like contributions) have destructive interference, and the iJ-resonance and box diagrams 
have constructive (destructive) interference for Mhh < Mh {Mhh > Mh)- In the SM, 
the h-resonance and box diagrams have exact destructive interference at the double-Higgs 
threshold. In the singlet model the cancellation is not exact anymore due to changes in 
the tri-linear coupling and different mixing angle suppressions of the two diagrams, but 
the SM like contributions still have the strongest destructive interference at threshold. For 
Mh < 2m/i, both h- and //-resonance diagrams have strong destructive interference with the 
box diagram near Mhh ~ 2mh. Hence, the overall destructive interference dip is strongest 
for Mh = 200 GeV. 

In picking a parameter point, we have a choice as to whether to choose a positive or 
negative sign for sin 6*. The comparison of these two choices is shown in Fig. [HI with the 
LHS showing the differential cross sections and the RHS the ratio of total cross sections. 
As shown in the LHS of Fig. [HI the choice of sign makes little difference in the shape of 
the distributions. In particular, the interference effects remain essentially unchanged. This 
can be understood by analyzing the triple couplings Am (Eq. [10]) and A 211 (Eq. [TT]), and 
Fi (Eq. [ 20 I) . The dependence of the cross section on the sign of sin/' always appears with 
an associated factor of tan [3 and is suppressed compared to the cos 9 terms in the triple 
couplings. However, there can still be a signihcant change in the total rate, as shown on 
the RHS of Fig. [HI For Mh > ^ruh, the cross section for negative sin/' is ~ 70 — 80% that 
of the cross section for positive sin/'. For Mh < ^.mh, the two cross sections are nearly the 
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pp—>hh (Singlet Model), VS = 13 TeV pp—>hh (Singlet Model), Vs=13 TeV 

tan p = 1.0, cos 9 = 0.9, |i = CT12NLO PDFs tanp=0.5. cose=0.96, CT12NLO PDFs 



FIG. 7: Exact LO rates for pp —)• hh at ^fS = 13 TeV for fixed singlet mixing parameters, 
cos 6* = 0.9 and tan/3 = 1.0 (LHS) and cos 9 = 0.96 and tan/3 = 0.5 (RHS). 

same. This can be understood, and is shown later, by noting that the H resonance makes 
a subleading contribution for Mh < 2m/j GeV and the SM like contributions only depend 
on the sign of sin0 in a highly suppressed sin^^ term in Am. Throughout the rest of the 
paper we will choose a positive sign for sin0. 

In Fig. [9l we show the ratio of the singlet model rate normalized to the SM rate. It is 
clear that near the resonances large enhancements in the rates are possible and the singlet 
model should be clearly distinguishable from the SM. 


C. Interference effects 


The presence of the second scalar leads to interesting interference effects with the SM-like 
contributions. The real parts of the propagators in (Eq. [20|) interfere destructively for 
rrih < Mfih < Mh and constructively for > Mh, as is typical for the interference of two 
resonances^. However, in the SM the box and triangle diagrams destructively interfere, with 
the box diagram dominating at large Hence, although the propagators of the two 

resonances destructively interfere, the TT-propagator constructively interferes with the box 
diagram for M^h < Mh, and destructively interferes for M^h > Mh- 

Leading order differential cross sections with individual contributions are shown sepa- 


This same interference effect is seen in the process gg ^ ZZ hv the singlet model 53|, and in Drell-Yan 


production below the Z peak. 
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pp^hh (Singlet Model), Vs=13 TeV pp^hh, Vs = 13 TeV 



FIG. 8: Exact LO rates for pp —)• hh at \fS = 13 TeV for fixed singlet mixing parameters 
cos 6 = 0.96 and tan/3 = 0.5. (LHS) The solid (dashed) curves correspond to choosing a positive 
(negative) sign for sin0. (RHS) The ratio of the total LO cross sections with negative and positive 
sign for sin0. 


pp^hh (Singlet Model/SM), Vs=13 TeV 



FIG. 9; Exact LO rates for pp —)• hh normalized to the exact LO SM rate at y/S = 13 TeV for 
fixed singlet mixing parameters, cos 6 = 0.96 and tan/3 = 0.5. 

rately in Figs. [10] and [TT] The curves labelled “h + H resonances only” include the con¬ 
tributions of both s-channel h and H and their interference, but not the effect of the box 
diagrams. The destructive interference between the two propagators for Mhh < Mh is clear. 
The curves labelled “no TT-resonance” have the H resonance contribution removed; that is, 
only the SM-like contributions are included. As described above, by comparing the curves 
labelled “no H-resonance” with the total distribution, we see that there is constructive in- 


14 

















pp^hh, VS= 13 TeV pp^.hh, Vs=13 TeV 

Mj^=200 GeV, CT12 NLO PDFs, cose=0.96, tanp=0.5 Mj^=300 GeV, CT12 NLO PDFs, cos0=O.96, tanp=0.5 



FIG. 10: LO results for pp ^ hh at VS = 13 TeV for Mh = 200 GeV (LHS) and Mh = 300 GeV 
(RHS), cos 9 = 0.96 and tan/3 = 0.5. See text for description of individual curves. 

pp—>hh, VS=13 TeV pp^hh, 3s=13 TeV 

Mj^=400 GeV. CT12 NLO PDFs, co.s0=O.96, tanp=0.5 Mj^=600 GeV, CT12 NLO PDFs, cos0=O.96, tanp=0.5 




FIG. 11: LO results for pp ^ hh at VS = 13 TeV for Mh = 400 GeV (LHS), Mh = 600 GeV 
(RHS), with COS0 = 0.96 and tan/3 = 0.5. See text for description of individual curves. 

terference between the H and SM-like diagrams for Mhh < Mh and destructive interference 
for Mhh > Mh- It is apparent that the —?■ oo limit fails to reproduce the correct interfer¬ 

ence structure near and slightly above the peak and overshoots the rate at high Mhh- The 
location of the interference dip just above the resonance is slightly shifted to larger Mhh in 
the —)■ oo limit. This motivates weighting the NLO rate (which is only known in the 

—)■ cx) limit), by the exact LO rate. 

We show the ratio of the interference between the H resonance and SM-like diagrams and 
the full invariant mass distribution in Fig. [121 Exact rrit dependence has been kept. The 
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pp^hh (Singlet Model), Vs = 13 TeV pp^hh (Singlet Model), VS = 13 TeV 

tan p = 1, cos 9 = 0.9, |X = CT12NLO PDFs tan P = 0.5, cos 0 = 0.96, CT12NLO PDFs 



FIG. 12: Ratio of the interference between the Fl-resonance and the SM-like contributions, h- 
resonance+box, and the full invariant mass distribution with tan/3 = 1.0 and cos 9 = 0.9 (LHS), 
and tan/3 = 0.5 and cos 9 = 0.96 (RHS). The curves use the exact LO matrix elements. 


interference contribution is 

d(^lo ^ daio _ f da^o 

dMhh dMhh \dMj,n dMhh )' ^ ^ 

where contains only the contribution from the i^-resonance, and contains the 

h-resonance and box contributions and their interference. An interesting feature of Fig. [12] 
is that for M^h ^ the interference contribution is independent of for hxed 6 and 
tan (3. This somewhat surprising effect can be understood by taking Fi (Eq. |20|) in the limit 
ml,s ^ Mjj: 


Fi ^ 


cos OXiiiv 
- ml + imhTh 


+ 


sin 6 sin 26 


cos6* + sin0tan/3) 


+s cos^ 9Fci{s, t, u, ml). 


(34) 


As can be clearly seen, in this limit, the double Higgs rate does not explicitly depend on the 
heavy scalar mass. 

The ratio of the interference between the iJ-resonance and SM-like contributions dehned 
in Eq. [33] and the total cross section are shown in the LHS of Fig. [13] We also show 
the ratio of the Ff-resonance contribution only and the total cross section in the RHS of 
Fig. [I3] The curves are shown for the two parameter points tan/3 = 1, cos 6 = 0.9 (solid 
black) and tan/3 = 0.5, cos9 = 0.96 (dotted red). At amplitude level, the dominant (box) 
contribution to the SM-like pieces is proportional to cos^ 9 and makes a similar contribution 
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for both parameter points. However, below 2mh the hT-resonance amplitude is proportional 
to sin^ 9 and sensitive to relatively small changes in cos 9. This explains why for Mh < 
2mh the interference and TT-resonance contributions are larger for cos6' = 0.9 than for 
cos^ = 0.96. For > 2mh and using the narrow-width-approximation, the iJ-resonance 
amplitude is proportional to sin0 and is still larger for cos6' = 0.9 than for cos 9 = 0.96. 
Once the resonance production of hh turns on, Mh ~ 2mh, the FT-resonance contribution 
dominates, as seen in the RHS of Fig. [131 As Mh increases, the iF-propagator suppresses 
the Ff-resonance contribution. However, as Mh approaches 2mt, as is well-known in single 
Higgs production, the production rate through a top quark triangle increases. For 2mh ^ 
Mh ^ 2mt these two effects cancel each other and the contribution from the Ff-resonance is 
relatively constant. As Mh increases above ~ 2mt, the suppression from the Ff-propagator is 
the dominant effect. Hence, the fractional contribution from the Ff-resonance only decreases 
and the fractional contribution from interference increases. These two effects are correlated 
because the SM-like contribution by itself is independent of Mh- It should be noted that the 
absolute contribution from the interference is nearly independent of Mh for Mh ^ 500 GeV. 
This can be understood from Eq. |33l Since for increasing Mh there is a large contribution 
to the cross section from the Mh ^ M^^ region, the total contribution to the interference 
is largely independent of Mh- 

D. NLO Effects 

In Fig. im we show the enhancement of the total cross section in the singlet model, 
relative to the SM rate. For tan/3 = 5 and cos9 = 0.96, the maximum enhancement is of 
(9(8) for Mh ^ 500 GeV and decreases rapidly to (9(1) for larger Mh- For larger mixing, 
tan /3 = 1 and cos 6* = 0.9, enhancements of the SM rate up to a factor of ~ 22 are possible. 
We see that jasM is not very different for LO and NLO total rates. The contribution 
of the H resonance in the narrow width approximation is accurate for Mh ^ 400 GeV^ but 
underestimates the enhancement for larger Mh- 

We now present our numerical results for the double Higgs invariant mass distributions at 
NLO. Fig. IT5] shows the individual contributions (Eq. l25ll to the invariant mass distributions 
using the approximation of Eqs. 12511261 It is important to remember that the full rrit 
dependent NLO rate is not known. We plot the absolute value of the qg contribution, since 
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pp^hh (Singlet Model), Vs = 13 TeV pp^hh (Singlet Model), VS = 13 TeV 

ji = CT12NLO PDFs |i = CT12NLO PDFs 



FIG. 13: Ratio of the interference between the ff-resonance and the SM-like contributions, h- 
resonance+box, and the total cross section as a function of Mh (LHS). Ratio of the H-resonance 
contribution only and the total rate (RHS). Both tan/3 = 1.0, cos 9 = 0.9 (solid black) and tan/3 = 
0.5, cos 9 = 0.96 (dotted red) are shown. The curves use the exact LO matrix elements. 


pp^hh, Vs = 13 TeV 

^ = CT12NLO PDFs 



FIG. 14: Ratio of NLO (dash-dot black) and LO (dashed red) rates to the respective SM rates. 
Also, the ratio of the LO H-resonance only contribution calculated using the Breit-Wigner resonance 
(solid green) and in the narrow-width approximation (blue dotted) to the LO SM rate. 

it is negative. The leading corrections are from the gg and virtual contributions, while the 
qg and qq contributions are subleading. 

It is interesting to compare the effect of the approximations to the top mass dependence 
at NLO. In Fig. [161 we compare the NLO rate for Mh = 300 GeV computed using the 
approximation of Eqs. I2S1I2H] (dashed red curve) with that obtained by computing 
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pp^hh (Singlet Model), Vs=13 TeV 

M^=300 GeV, taiiP=0.5. cose=0.96, CT12NLO PDFs 



FIG. 15: Total and individual contributions to the NLO cross sections defined in Eqs. 1251261 The 
model parameters were chosen to be tan/3 = 0.5, cos0 = 0.96, and Mh = 300 GeV. 

(Eq. [30]) and weighting by the exact rrit dependent LO cross section (solid black). The curves 
overlap almost exactly. Since most contributions to the NLO rate (Eq. [25]) are proportional 
to the LO rate, the approximate rrit dependence is mostly captured by weighting the exact 
LO rate with The only complication is a piece of the virtual contribution (Eq. |26|) 

that is not proportional to the LO rate. However, this piece turns out to make a small 
contribution. 

We then compare with an NLO rate computed in the rrit ^ oo limit (dotted blue in 
Fig. [T6|) . i.e. this result is not reweighted by the exact rrit dependent LO result. The 
rrit ^ oo limit shifts the location of the interference dip to slightly higher Mh^. This effect 
is also apparent in the comparison of the exact rrit dependent and rrit ^ oo LO curves of 
Fig. [TOl A blow-up of the interference region is shown on the RHS of Fig. [I6] and makes 
this effect obvious. 

On the RHS of Fig. [161 we can also see that the curve obtained by weighting the exact 
LO rate by differs from the curve calculated using Eqs. 125112^ at the interference dip 

of the rrit ^ oo curve. The interference dip is where the LO cross section is a minimum. 
Hence, the piece of the virtual contribution (Eqs. I261I3T]) not proportional the LO cross 
section makes a relatively large contribution in this region. Since the interference dip is 
deeper in the rrit ^ oo limit (see Fig. ITOll . this effect is more pronounced in the rrit ^ oo 
case. As a consequence, at the interference dip, the rrit ^ oo NLO rate is not approximately 
proportional to the rrit ^ oo LO rate. Therefore, weighting the exact LO rate by 
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pp-4hh (Singlet Model), VSj,=13 TeV pp^.hh (Singlet Model), VSjj=13 TeV 



FIG. 16; NLO cross sections for Mh = 300 GeV and with different approximations for the top 


mass dependence as described in the text. 
cos 9 = 0.96. 


pp—>hh (Singlet Model), VSjj=13 TeV 



mixing parameters were set to tan /3 = 0.5 and 


pp^hh (Singlet Model), a/Sjj= 13 TeV 



FIG. 17: NLO cross sections for Mh = 200 GeV and with different approximations for the top 
mass dependence as described in the text. The mixing parameters were set to tan/3 = 0.5 and 
cos 9 = 0.96. 

does not reproduce the curves computed using Eqs. 12511261 precisely where the rrit ^ oo rate 
has the strongest destructive interference. 

It is interesting to compare with the NLO rate for a heavy Higgs mass below the threshold 
for a double Higgs resonance, Mh = 200 GeV. These results are shown in Fig. [171 In the 
interference region, the effects are similar, but more pronounced, to those in the Mh = 
300 GeV case. In fact, for Mh = 200 GeV, the two curves computed by Eqs. I25|26l and by 
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FIG. 18: Ratio of differential K-factors evaluated in the nrit ^ oo limit to those calculated using 
the approximate mt dependence of Eq. [25l The mixing parameters were set to tan (3 = 0.5 and 
COS0 = 0.96. 


weighting the exact LO rate by do not agree at the minimum of the curve in 

addition to the minimum of the curve. This can be understood by noting that 

as Mh increases, the interference dip of the LO cross section is more shallow (see Fig. [7j). 
As a consequence and discussed above, as Mh increases the contribution to a^rt that is 
not proportional to the leading order rate decreases. Hence, the curves computed using 
Eqs. 12511^ and weighting the exact LO rate with will be in better agreement with 

increasing Mh- In Fig. [18] we show the ratio of the and (Eq. [29]), which is the 

same as the ratio of the NLO rates calculated by weighting of the exact LO rate by 
and using Eqs. I25][26l As can be seen, as Mh increases the two methods increasingly agree. 

In Fig. [ini we show the scale dependence of the invariant mass distribution for a represen¬ 
tative parameter point with Mhh/2 < < 2Mhh. The LO cross sections contains exact rrit 

dependence and the NLO cross section is computed using Eqs. [25][26l The NLO corrections 
decrease the scale dependence from ~^ 2 o% ~ ±15%. Additionally, the NLO scale 

dependence is fairly flat throughout the distributiou; in particular, it does uot appreciably 
chauge iu the resouauce aud stroug destructive interference regions. 

In Fig. [20] we show the differential K-factor in the —)■ cxo limit, as dehned 

in Eq. [30] The K-factor is flat with a value of 2 — 2.2, except for spikes that occur in 
the regions with the strongest destructive interference. As shown in Fig. [TH] the K-factor 
computed using Eq. [29] agrees with except in the regions of strong destructive 


21 


















pp—>hh (Singlet Model), VS = 13 TeV 

tan p = 0.5, cos 0 = 0.96, M^^/2 < 0 < 2 CT12NLO PDFs 



pp—>hh (Singlet Model), VS = 13 TeV 

tan P = 0.5, cos 0 = 0.96, < p < 2 CT12NLO PDFs 
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FIG. 19: Scale dependence of the invariant mass distributions for Mh = 300 GeV for both the NLO 
and LO cross sections with Mhh/‘^ < ^ < 2Mhh- The RHS shows the fractional scale dependence of 
the invariant mass distributions. These hgures are computed using Eq. [25j The mixing parameters 
were set to tan/3 = 0.5 and cos 6 = 0.96. 


pp^hh (Singlet Model), VS = 13 TeV 

tan p=0.5, CO.S 6=0.96, CT12NLO PDFs 


pp—>hh (Singlet Model), VS = 13 TeV 

tan p=0.5, cos 6=0.96, CT12NLO PDFs 




FIG. 20: as defined in Eq.[29]at y/S = 13 TeE for various Mh and tan (3 = 0.5, cos 9 = 0.96 


interference. 


E. Results at 100 TeV 

Next we present onr resnlts for the NLO calculation of double Higgs production at y/S = 
100 TeV. In Fig. [2l] we plot (top) Eq. |30l and (bottom) the ratio of the 

and iF™*, Eq. [29l The K-factors at y/S = 100 TeV are similar to those at y/S = 13 TeV. 
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pp->hh (Singlet Model), Vs = 100 TeV pp^hh (Singlet Model), VS = 100 TeV 



M„„(GeV) M^JGeV) 


FIG. 21: (Top) defined in Eq. [30l (Bottom) Ratio of and R'™'* defined in Eq. [29l 

The mixing parameters were set to tan/3 = 0.5 and cos0 = 0.96 

Since the ratio of K-factors at 100 TeV is similar to those at 13 TeV, our comparison of the 
rates calculated by weighting the exact LO rate by and Eqs. will translate 

from the 13 TeV to 100 TeV environment. 

In Fig. [22] we show the normalized invariant mass distributions at y/S = 13 and 100 
TeV with both —)■ cxo and the approximate hnite rrit dependence of Eqs. 12511261 As 
noted previously, the inhnite top quark mass limit overestimates the tail of the distribution. 
Additionally, for the SM-like contributions, the —)■ cxo limit underestimates the cross 

section for M^h ^ 550 GeV (Fig. 0]). Hence, after the strongest destructive interference, 
the SM-like contribution to the approximate finite rate grows more quickly than in the 
—)■ cx) case. As a result, directly after the interference dip, the approximately finite mt 
distribution grows more quickly and obtains a higher value than the f oo distribution. 


23 


























pp—>hh (Singlet Model) 

= 300 GeV, lan p = 0.5, cos 0 = 0.96, Jl = M^|_. CT12NLO PDFs 



FIG. 22: Normalized invariant mass distributions for ^/S = 13 and 100 TeV in both the —?> oo 
limit and the approximate mt dependence of Eal25l The model parameters were set to Mh = 
300 GeV^ tan/3 = 0.5 and cos 9 = 0.96. 


Finally, at \fS = 100 TeV the tails of the distributions are enhanced relative to 13 TeV. 
This is because for a given invariant mass, the PDFs are evaluated at smaller x at 100 TeV 
than at 13 TeV. Hence, the enhancement of the gluon parton luminosity causes the tail of 
the distribution to be longer. 


V. CONCLUSIONS 

The production of Higgs pairs from gluon fusion is an important probe of the structure of 
the scalar potential. In the SM, the QCD corrections are known in an approximation where 
the LO rate is weighted by a iP factor computed in the —)■ cxo limit, increasing both the 
total rate and da/dM^h by a factor of around 2. 

We have presented results in the Higgs singlet model, where the tri-linear Higgs self 
coupling is modihed from the SM value and significant resonant effects from the second scalar 
occur. The effects of the interference between the heavy scalar and SM-like contributions 
can be significant, altering invariant mass distributions for all M^- For Mh ^ 450 GeV, the 
interference effects can make a ~ 10—20% contribution to the total rate. For Mh ^ 2m/i, the 
interference effects can suppress the total cross section up to ~ 30% for a viable parameter 
point. Hence, in searches for heavy scalars, these effects should be included. 

We compare an approximation for the NLO QCD corrections where the exact mt de¬ 
pendent LO cross section is weighted by a iP factor computed in the —)■ cxo limit, and 
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alternatively where the exact rrit dependent form factors are inserted into the NLO contri¬ 
butions. The approaches give similar results except in the regions with large destructive 
interference. 

In the singlet model, the total cross section is increased by factors between 5 — 10 above 
the SM rate for tan/3 = 0.5 and cos6 = 0.96. For larger mixing (tan/3 = 1 and cos 6 ' = 0.9), 
we hnd enhancements from the SM rate between 10 — 20 for Mh < 500 GeV, and the 
enhancement is very similar at LO and NLO. The resonant approximation to the total cross 
section underestimates the enhancement by about a factor of 2 at large Mh- 

The singlet model demonstrates a case where the kinematic distributions of the outgoing 
SM Higgs pair are signihcantly altered from the SM, and where the higher order QCD 
corrections differ from those of the SM near the resonance peak. 
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